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We discuss the recent theory of Gvozdikov [Phys. Rev. B 70, 085113 (2004)] which aims at 
explaining the Shubnikov-de Haas oscillations of the longitudinal resistivity pzz observed in the 
quasi-two-dimensional organic compound — (BEDT — TTF)2SF5CH2CF2S03. We point out that 
the self-consistent equations of the theory yielding the longitudinal resistivity and the magnetic field 
dependence of the chemical potential have been incorrectly solved. We show that the consideration 
of the self-consistent Born approximation (which determines the relaxation rate in Gvozdikov's 
paper) leads in fact to the complete absence of the longitudinal conductivity azz at leading order in 
high magnetic fields. 

PACS numbers: 72.15.Gd, 75.20.En, 73.43.Cd, 71.18.-l-y 
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Gvozdiko\»i recently derived a theory of the Shubni- 
kov-de Haas oscillations of the longitudinal conductivity 
azz in quasi-two-dimensional (2D) metals in high perpen- 
dicular magnetic fields. He found that azz minima dis- 
play a thermally activated behavior and that azz peaks 
are split in presence of small oscillations of the chem- 
ical potential /i. He also claimed that his theory can 
explain the giant oscillations of the magnetoresistance 
Pzz = a~^ observed^ in the layered organic compound 
P" - (BEDT - TTF)2SF5CH2CF2S03. 

In this Comment, we point out that the theoretical 
results have been incorrectly derived and that the theo- 
retical model adopted by Gvozdikov is inconsistent with 
existing experiments. In the following, we first address 
the framework of the theory and the problems with the 
calculation of azz- Then, we shall focus on the thermo- 
dynamic quantities such as the density of states and the 
chemical potential. Finally, we shall discuss possible the- 
oretical explanations of experimental observations. 



I. THE LONGITUDINAL CONDUCTIVITY a^z 
A. Model and framework 

The system under consideration is a quasi-two- 
dimensional metal with an interlayer hopping integral t. 
In the paper the interlayer dispersion is modelled by an 
unknown density of states with a characteristic width of 
the order of t. Gvozdikov's paper— was aimed at deriv- 
ing the oscillations of azz as a function of magnetic field 
perpendicular to the layers under the regime of a weakly 
incoherent interlayer transport. 

The incoherent transport was discussed in literature 
(see, e.g., Ref. 3) in the cases of the so-called weakly and 
strongly incoherent regimes. The former occurs when 
the intralayer momentum is conserved in the process of 



transfer of the electron from one layer to another. The 
incoherence is relevant when a large number of in-plane 
scatterings takes place before an interlayer tunneling oc- 
curs, i.e. when t <^ fi/T, where r is the mean free time. 
Interlayer transport is then incoherent because the suc- 
cessive tunneling events get uncorrelated. This situation 
was first discussed in the paper* in relationship with the 
conductivity in cuprates and it gives a Drude-type for- 
mula for the interlayer conductivity where the scattering 
time is the time of inelastic intralayer scattering. The 
strongly incoherent transport occurs when the intralayer 
momentum is not conserved by tunneling and there is 
no interference between the wavefunctions on adjacent 
layers. 

In Gvozdikov's theorjsi intralayer inelastic processes 
are apparently^ not considered since the width F = h/2T 
of the spectral function in the layer is determined by the 
elastic scattering on impurities (note that in our defini- 
tion there is a difference of a factor 2 with Gvozdikov's 
notation). The additional complicacy considered in Ref. 
[l| consists of the inelastic processes at tunneling. The 
Landau levels and (implicitly) the centers of the Landau 
orbits are assumed to be conserved. This is valid pro- 
vided that the characteristic value of the energy change 
at tunneling, 7, is much less than the intralayer spectral 
width F. Under this regime of so-called weak incoherence 
by Gvozdivov (that does not correspond to the definition 
of weak incoherence given above) , the inelastic processes 
at tunneling do not play a role for the interlayer transport 
and we return to a theory where only elastic scattering 
on impurities is relevant. It should be noted that there 
is no clear distinction between t and 7 in Ref. [l|. 

The derivation of the longitudinal conductivity azz for 
elastic impurity scattering and a cosinclike interlayer en- 
ergy dispersion was already reported in the Rcfs 6.7,^. 
The paper— was interested in the regime hfl < t (here 
il is the cyclotron frequency), while the papers^i^ were 
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concerned with the 2D limit Ml t in connection with 
the experiment)^ 



B. Derivation of a^z 

The longitudinal conductivity a^z is derived within the 
Kubo formula written in the formalism of the Green func- 
tions. Vertex corrections have been neglected in Ref. 1 
which means that a point-like impurity model has been 
implicitly assumed. Within this model the self-energy 
part entering into the expression of the disorder averaged 
Green functions is independent of the Landau levels in- 
dex n. This fact allows to compute analytically the sum 
over n using the Poisson summation formula. 

Gvozdikov has derived an expression for azz, Eqs. 
(21)-(22). In addition to the condition t < T, the 2D 
limit h^l 3> i is required for its validity. Under the same 
inequalities, one can rewrite these equations (21)-(22) as 



I dE{-f'{E))[aB 



where 
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Here f{E) is the Fermi-Dirac distribution function which 
depends on the chemical potential /i, ar is the conductiv- 
ity at zero magnetic field, X{E) = n/i}T{E), Xq = tt/VLtq 
[t(E) and tq are respectively the elastic quasi-particle 
lifetimes at energy E in the presence and absence of mag- 
netic field] and the function S is defined by 



S{X{E),E) 



smhX{E) 



coshA(£') +cos(27r£'/rifi)' 



(4) 



Formulas (l)-(4) are incomplete without a prescription 
for calculating the function X{E). To proceed further, 
Gvozdikov applied the self-consistent Born approxima- 
tion (SCBA) to the quasi-2D spectrum. This approxima- 
tion stipulates the proportionality between the density of 
states N{E) and t-'^{E), i.e., 
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Here A'o is the density of states of the quasi-2D system 
at zero magnetic field. The density of states N{E) under 
magnetic field is expressed as a Fourier series and is a 
function of the lifetime t{E). In the 2D limit {hil ^ t) 
and for t <C F, it is possible to sum up this series and the 
relation (5) leads to the following self-consistent equation 
for X{E): 



X{E) = XoS{X{E),E). 



Equation ([6]) yields straightforwardly the fact that the 
term as = cr^ [see Eq. i.e., this latter does not os- 
cillate with the magnetic field. Gvozdikov got then that 
the oscillations in the strong 2D regime arise from the 
term ctq only. Here we want to stress that this result- 
has been incorrectly derived from the self-consistent Eq. 
([S]). All the subsequent results based on the formula (28) 
for azz in Ref. [l| are accordingly incorrect. Indeed, dif- 
ferentiating the formula ([6]) with respect to A, we obtain 
that XodS/dX — 1. Substituting this result into Eq. ^ 
we straightforwardly have ctq = —a^, which means that 
in frame of SCBA the total conductivity 







(7) 



in the limit HQ 3> t. A non-zero expression for azz may 
be obtained by considering next order corrections in the 
parameters t/Hfl and t/T. 

It is worth noting that the Eqs. are exactly the 

same as Eqs. (19)-(21) given in the work! considering 
a cosinelike dispersion relation with respect to the inter- 
layer momentum. The reason for this is that the specific 
form of the interlayer density of states is unimportant 
under the conditions t <C hfl and t <^ H/t. The impor- 
tant difference between the two theories is that in Ref. 
the total (not just the intralayer) spectral width due to 
elastic scattering on impurities F = h/2T is postulated 
as being constant, i.e. energy independent: the SCBA 
leading to Eq. ([6]) and Eq. ^ is not considered as an 
inherent property of the transport theory. 

We would like to stress that Eq. ([7]) is a consequence of 
the application of SCBA in high magnetic fields only. In 
smaller magnetic fields < t, the same SCBA yields^ 
small oscillations of azz around the zero field value which 
is nonzero due to 3D elastic scattering on impurities. 



II. THERMODYNAMIC QUANTITIES 
A. Density of states 

It is important to note that Eq. ^ determines as well 
the function X{E) (or the hfetime t{E) — n/flX{E)) as 
the density of states N{E) — NoX{E)/Xo. It has been 
recognized in the Section IV of Ref. that the self- 
consistent equation ^ for the density of states is ex- 
actly the same equation as one would obtain using the 
SCBA with a strict 2D spectrum, which has has been 
well known for a long time.— The approximation con- 
sisting in replacing X{E) everywhere by the peak value 
X{En) ~ V2Ao as done in Ref. [l| is inconsistent with Eq. 
(I6|). In fact, the self-consistent Eq. ([6|) yields a function 
X{E) which strongly oscillates with E, especially in high 
magnetic fields. When Aq < 2, X{E) even vanishes be- 
tween Landau levels energy positions £"„ = {n + \/2)Ml, 
which means that the density of states splits int o sep arate 
bands centered on the energies En (see Refs. iEl). Go- 
(6) ing beyond SCBA, it is possible to calculate more accu- 
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rately the tails of the bands to avoid the unphysical cut- 
off produced by Eq. Moreover, next corrections due 
to the interlayer hopping, of the order of t/hil. and t/T, 
may yield finite but small values for N{E) between Lan- 
dau levels. However, in any cases, the model adopted by 
Gvozdikov automatically implies sharp peaks for N{E) 
when E w and the quasiabsence of states between 
Landau levels for Qtq > tt/2. 



B. Chemical potential 

Eq. (23) of Ref. 1' [reproduced here with Eq. ([6])] com- 
pletely determines the density of states and thus allows 
to compute the chemical potential /i as a function of mag- 
netic field. Therefore, there is no freedom for the choice 
of some phenomenological model describing the chemical 
potential oscillations in accordance with the experiments, 
as surprisingly done in the Section 4 of Ref. [l|. 

If the density of states consists of sharp bands centered 
on the Landau levels, one expects that fi is pinned to a 
value close to (n -|- l/2)h^ for most of the range in mag- 
netic fields and drops suddenly to another Landau band 
once a Landau band becomes filled or empty. This pro- 
cess gives rise to strong magnetic quantum oscillations 
of the chemical potential. This expectation is physically 
inconsistent with the experimental observation of an in- 
verse sawtooth profile of the magnetization oscillations in 
the compound /?"- (BEDT - TTF)2SF5CH2CF2S03 ,— 
that impliesi^ the presence of negligibly small oscillations 
of fi. Furthermore, the experiments on thermodynamic 
and transport quantities^iii rather suggest a significant 
amount of states between the Landau levels for ^Itq ~ 1 
in contradiction with the result obtained within the self- 
consistent Born approximation. 

Now, we discuss the different forms possible for the 
chemical potential oscillations under the assumption of a 
constant lifetime (as in Section 4 of Ref. [ll), because they 
have been incorrectly presented and used in Ref. [J. For 
an energy independent lifetime f, the equation obeyed at 
zero temperature by the chemical potential fi in the 2D 
limit can be explicitly calculate d^^'^'^ 



n = Ep H arctan r— — 



(8) 



where v = Tr/fif and Ep is the Fermi energy at zero mag- 
netic field. As pointed out by Grigoriev,^'* it is possible 
to invert this self-consistent formula (O to obtain 



Ml ( sm{2TTEF/hn) 
Ll — Ep H arctan — — — 



(9) 



When <C 1, the density of states within this model 
consists of sharp peaks centered around the Landau lev- 
els (as in the SCBA), and this Eq. ^ yields a staircase 
dependence for fi/hfl as a function of magnetic field. Cor- 
respondingly, the oscillating part jl/Ml = {fi — Ep)/Ml 



exhibits a sawtooth dependence (with the so-called direct 
shape) 

Within the phenomenological model assuming a reser- 
voir of states and a constant lifetime, the equation for /i 
now takes the fornti^ 



^ f sm(27ru/ni]) \ , , 

Ep + — — arctan — ^— (10) 

tt{1 + R) ye" + cofi{2TTfM Ml) J ^ ' 



where i? is a dimensionless parameter measuring the 
strength of the reservoir. Equation pH]) embraces nat- 
urally the former model with i? = 0. For an important 
reservoir, i? ^ 1, the oscillating part \i is extremely small 
(/i <C Mt)^ implying that [i is almost fixed to the zero field 
value Ep. In this particular case, it is then justified to 
replace /i by Ep in the right-hand side of this Eq. ([TT 
to obtain 



Ml ( s\n{2TTEp/Ml) 

a w Ep H arctan r-— — 

^ TTi? \e'' +cos{2TTEp/hVl) 



(11) 



Then, jl/Ml exhibits the (so-called) inverse sawtooth 
waveform when <C 1, but with a negligibly small am- 
plitude (because i? ^ 1). On the contrary, the magneti- 
zation oscillations which are proportional to Rjl/Ml are 
not limited in amplitude by the reservoir parameter R 
(see Ref. [H). 

In his discussion of the effects of the chemical potential 
oscillations on cr^^, Gvozdikov^ has used Eq. ([5]) (i.e. no 
reservoir) with fi replaced by Ep in the right-hand side 
of the equation. We want to stress that it is generally 



mathematically incorrect to replace /i by Ei 
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the 



right-hand side of Eq. ([8]) especially when v <^ \ [e.g., 
compare the resulting equation and Eq. An addi- 

tional remark is that Eq. (31) for fi considered in Ref. 
m [Eq. with a minus instead of a plus between the 
two terms in the right-hand side] to describe the direct 
sawtooth case has no basis at all. 



III. THEORETICAL APPROACHES 

As already mentioned above, the theory in Ref. 0, also 
aimed at explaining the oscillations of azz reported in 
Ref. 0, is based on the same Eqs. (l)-(4) for Ozz- The 
difference with Gvozdikov's theory rests on the function 
\{E) [or equivalently on the lifetime t{E)], which is (phe- 
nomenologically) assumed to be constant. Within this 
model both contributions ub and ctq oscillate with mag- 
netic fields. The thermal activation of Uzz minima stems 
from a cancellation at zero order in t/Ml of these two 
contributions between the Landau levelsi^ Contrary to 
SGBA, the cancellation does not occur at the Landau 
levels. 

Now, we would like to address one alternative expla- 
nation to the experiment.^ Owing to some similitudes, 
it is tempting to bring together the present problem in 
a layered conductor and the problem encountered in the 
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theory of the quantum Hall eSec^. At the experimental 
side, in this quasi-2D conductor and in the 2D electron 
gases under the regime of the integer quantum Hall effect, 
a thermal activation is observed: in the first system it 
concerns the minima of the interlaycr conductivity azz r 
and in the second system the minima of the longitudinal 
conductivity (Txx-^^ The puzzling point is that the den- 
sity of states extracted from different measurements of 
the thermodynamic quantities points out a rather large 
amount of states between the Landau peaks in the lay- 
ered organic conductorSiii and in the 2D electron gases 
under the regime of the integer quantum Hall effec t^^i^^ 
which is typically observed when Qtq > 1. As a result, 
we have to face the same difhculty, i.e., to capture trans- 
port and thermodynamic properties within a complete 
and consistent microscopic theory. In the regime of the 
integer quantum Hall effect, it is believed that most of 
the states of the 2D layer are localized by the disorder 
and act as a reservoir which almost fixes the chemical 
potentialji^ The core of the extended states is located 
at the energies En- At finite but not too low tempera- 
tures, the main means of conduction when the chemical 
potential sits between the Landau levels is the thermal 
excitation of quasiparticles at the edges of the mobility 
gap.— We could presume the presence of similar localized 
states in the layers responsible for the thermal activation 
of azz minima at high magnetic fields in the compound 
(3" - (BEDT - TTF)2SF5CH2CF2S03. 

For the moment there exists no proper quantum cal- 
culations for axx in the quantum Hall effect regime;^ 
so that this scenario for azz only rests on pure quali- 
tative considerations. In fact, the principal theoretical 
difficulty is the microscopic treatment of the impurity 
scattering effects on the Landau levels quantization. For 
this purpose, it is worth noting that the consideration of 
a model of impurity potential with zero-range is unphysi- 



cal in high magnetic fields B since in the limit B ^ oo the 
magnetic length — [he/ eB)^/"^ always becomes smaller 
than the correlation radius of the potential. It has been 
shown^^ for the 2D systems that the SCBA is no more 
valid when the range of the potential exceeds Ib so that 
the investigation of impurity effects within the usual per- 
turbation theory becomes very complicated. However, it 
seems that even this ingredient which consists to consider 
a finite range for the impurity potential is not sufficient to 
explain the presence of a background reservoir of statesi^ 
Apparently, new ideas are needed. 



IV. CONCLUSION 

The results reported in the paper— are based on the in- 
consequent use of SCBA. The proper SCBA application 
at leading order in high magnetic field leads to the zero 
conductivity azz = and to strong oscillations of the den- 
sity of states and of the chemical potential. Both latter 
effects can be suppressed in the presence of a significant 
number of states between Landau levels. An explanation 
for the thermally activated behavior of azz minima could 
be the presence of localized states as in the integer quan- 
tum Hall effect. However, a concrete self-consistent mag- 
nctotransport theory going beyond SCBA and including 
these localized states is still missing. 
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